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Abstract 

Consider the nth iterated Brownian motion = Bn o ■ ■ ■ o Bi. Curien and Konstan- 
topoulos proved that for any distinct numbers ti ^ 0, converges in 

distribution to a limit I[k\ independent of the tds, exchangeable, and gave some elements 
on the limit occupation measure of /F). Here, we prove under some conditions, finite di¬ 
mensional distributions of nth iterated two-sided stable processes converge, and the same 
holds the reflected Brownian motions. We give a description of the law of I[k], of the finite 
dimensional distributions of , as well as those of the iterated reflected Brownian motion 
iterated ad libitum. 
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1 Introduction 

Let B, Bi, B 2 ,.. ■ be a family of i.i.d. independent two-sided Brownian motions (BM), meaning 
that for any n, {Bn{t),t > 0) and {Bn{—t),t > 0) are two independent standard linear BM. 
Denote by = Bn o ■ ■ ■ o Bi the nth time iterated BM. Curien and Konstantopoulos [7] 
obtained the following results, gather in the following proposition. 

Proposition 1.1. (1) For any k > 1, any non zero ti,... Uk, converges 

in distribution. The limit distribution yik does not depend on the U’s, and then is exchangeable. 

(2) For {Ii,..., Ik) ^ ^k, the equality {Ii,..., U) ,B{Ik)) holds. Moreover, 

{I 2 — h, ■ ■ ■ ,Ik — Ii) — (R, • • • , Jfc-i) ~ hk-i- 

The distribution of Ii possesses the density exp(—2|a:|) over M (this result appeared first in 
Turban JT^l. 

(3) Let fin be the occupation measure of on [0,1], then the sequence {fin,Ti > 0) converges 
as n —)> 00 in distribution to a random probability measure fi, which has a.s. a finite support, 
and which has a.s. a Holder continuous density with exponent 1/2 — e for all e > 0. 

In this paper we go on this study in several connected directions; among other we give 
some elements on /r^, study iterated reflected BM, iterated stable processes, and provide a 
description of the finite dimensional distribution of the nth iterated BM 
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Here are the main lines of the paper. In Section we present the studied processes and fix 
some notations. In Section we provide some common features of the processes we iterate. 
Given a finite set of points L = {ii,i = 0, ...,/c}, the gaps sequence of L is the sequence 
G = {ii — ii-i, 1 < i < k), oi differences of successive points in L when sorted in increasing 
order. It turns out that for processes X with independent and stationary increments, the 
distribution of the gaps sequence of X{L) = {X{li)^i G /} can be described uniquely using 
the gaps sequence G of L. This simple property will appear to be at the heart of our advances 
about iterated BM. 

In Section devoted to iterated BM and iterated reflected BM, it is explained that if 
the initial gaps sequence G is a /c tuple of independent exponential random variables (r.v.) 
with parameters (Ai,...,Afc) then the gaps sequence of X[L) is distributed according to a 
mixture of k-tuple of independent exponential r.v., whose parameters are explicit functions of 
(Ai,..., Afc). To encode this property, we define a Markov chain > 1) at the parameter 

level, which makes explicit this parameter evolution (see (31) and around). 

A consequence is that the gaps sequence of the iterated BM ad libitum is a mixture of 
independent exponential r.v., and this mixture can be descr ibed precisely using the invariant 
distribution of the Markov chain (Propositions 


4.5 


4.6 


and Theorem 


4.10). 


Somehow, Remark 4.9 implies that our description of the iterated BM finite dimensional 
distribution, while complex, is the simplest we could expect. 

The same construction, using an analogous of the parameter Markov chain (Z^”\n > 1), 
implies that the law of the nth iterated BM is accessible if the gaps sequence of the initial 
distribution follows some independent exponential r.v. In Section 4.3, it is seen that this 


property provides a Lapl ace t ype transform of the finite dimensional distributions of the nth 
iterated BM Section ^ is devoted to the iteration of reflected BM. Section [H] is devoted 
to the iteration of stable processes, whose study appear much similar to that of iterated BM, 
except that explicit computations are out of reach for the moment. 

We discuss in Section [6] some natural extensions of this work. 


2 Random processes 

“BM” will be used to denote the two-sided linear BM as defined at the beginning of Section 
The process corresponding to the nth iterated process will be denoted (the process iterated 
under discussion, denoted X further, will be clear from the context). The processes iterated 
ad libitum, the limit of in the sense of the topology of finite dimensional distribution 
convergence, when it exists will be denoted I. The reflected BM is the (one-sided) process 
(|i?(t)|,t > 0) where B is the standard linear BM. 

We go on discussing stable processes (see Applebaum [T] for more information). We will 
consider only two-sided stable variables Z that can be written under the form A -|- r where 
A is stable symmetric (null skewness), and r a real number (the location parameter). The 
characteristic function of such a r.v. Z can be written under the form 

v^(u) = E(e*“^) = 


where 


ri{u) = — |tt|“ (t“ -|- iru 
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where a G (0,2] is the index of stability, a G (0,oo) the scale parameter (Theorem 1.2.21 in 
[I]). A stable process > 0) with parameters (a, cr, r) is the process such that Ai(0) = 0, 

with stationary and independent increments, and whose characteristic function is 


=E 


^iuX (t) 




( 1 ) 


The two sided stable process G M) is the process such that {X{t),t > 0) and {X{t),t < 

0) are independent and {X{t),t > 0) and {—X(—t),t > 0) are both one-sided stable process 
with parameters (a, a, r). For any t G M*, 


Xt - tr (d) 

|t|V« “ 


Ai-r. 


( 2 ) 


For any c > 0, {X{c°‘t),t > 0) is a stable process with parameters (a, ca, c“r). 

Let {Xi,X 2 ,...,) be a family of i.i.d. two-sided stable processes with parameters (a, a, r). 
The nth iterated stable process of parameters (a, cr, r) is the process 

/W =X,o...o Ai. 


We keep the same notation as for the iterated two-sided BM for some reasons that will appear 
clear below. 


Remark 2.1. The BM is the stable process with parameters (2, l/-v/2, 0). Its Markov kernel 
is F{Bt+s e dy\Bs = x) = exp(-(y - x)"^/{2t))/\/^. 


Iteration of stable processes with parameter (a, 1,0) and (a, a, 0) can be directly compared 


as explained in Remark 5.2 


3 Iteration of processes: general considerations 

In this section, we discuss some common features of the processes we iterate in the paper. 

All along the section k is a, positive integer: the size of the finite dimensional distributions 
under inspection. 

Notations. We denote by [a, bj the ordered sequence [a, b] n Z. The permutation group 
of the set [a, 6]| is denoted 5[[a, 5]|. Sometimes, we will use the notation x[a : 6] instead of 
{xa, ■ ■ ■ ,Xb), and also t[a : b],g[a : 6],A[a : b], etc, accordingly. The simple notation x[k] will 
stand for x[l : k]. 

For any sequence £[0 : k] = (^o, • • ■ denote by (.^[0 : k]) = sort(.^[0 : k]) this sequence 
sorted in increasing order. For any i G [[1, /cj, set 

The gaps sequence of i[0 : k] is the sequence of distances between the elements of {Iq, ■ ■ ■ ,1’fc}. 
It is defined by 

gaps(£[0 : k]) = G [[l,fc]|^ . 

Last, for x[l : k\ a sequence, x[0 : k] is the sequence defined by 

xo = 0,Xi = xi-\ - \-Xi, for i G [[1,/cj. (3) 
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Iteration of processes. What follows is valid for processes X such that X(0) = 0 a.s., 
with independent and stationary increments, which distribution are absolutely continuous with 
respect to the Lebesgue measure on M such that for any t > s, 

X{t) - X{s)= X{t-s) (4) 

whatever are the signs of s and t. Notice that this implies —X(s) X(—s) (taking t = 0). 
These general setting are satisfied by BM, by symmetric two-sided stable processes, and more 
generally, by symmetric two-sided Levy processes such that for any t > 0, X{t) owns a density. 
Some modifications are needed for processes as the reflected BM which have stationary but 
dependent increments. This is discussed in Section [4. 2[ 

Denote by <ht(.) the density of the distribution of X{t). We then have 

^tiy) = for any (t, y) e M* X M. (5) 

3.1 The gaps sequence evolution 

Let (to = 0,ti,...,tfc) be some distinct real numbers. We start with the description of the 
distribution of {X{ti),i G [[0, fej). As usual, the description is easier if the U are sorted... 

Let T G 5[[0,A:]] such that (L = ^ [0, fcj) = sort (t[0 :/c]). Hence , tr-i(o) = 0- 

Further let g[k] = gaps(t[0 : k]). The r.v. (A(tj+i) — X{ti),i G [I, A:]]) are independent, and 

X{ti) — X{ti-i) *"= X{Ati) depends on the gaps sequence of the L’s. Using the independence 
of the increments of X and their stationary, we get that the density / of {X{ti),i G [I, A;]]) on 
is 

k 

fiy[k]) = n i^yrU)) ■ (6) 

i=i 

where in the right hand side yo = 0. Indeed, one has {X{ti),i G [[0, kj) = {X{t^-Hi'^),i G [0, A;]]), 
and computing P (A(t.,—G dyi,i G [1, A:]]) = P {X{ii) G dyT-{i),i G [I, fcj) gives the result, 
using 0. 

The distribution of gaps(A(tj),i G [O, A:]]) depends also only on gaps(t[0 : A:]), and this is 
one of the key point of the paper. First, determine the vectors (A(tj),z G [O, A:]]) such that 

gaps{X{ti),iGlO,kj) = x[k] (7) 

some hxed element of (0, -|-oo)^. Clearly Q holds iff there exists some a G M such that 

sort(A(U),i G l0,k'l) = {a + Xi,i G [0, A:]]). (8) 

Equation Q implies that for a certain permutation r G 5|{0, A:]] 

iX{ti),i G [0,A:]]) = (a-hx^(i),i G l0,kj) 

from what we find 


{X{ti) - X{ti-i),i G [I, A;]]) = (Ax^(i),i G [1, A;]]). (9) 

The following proposition should be clear now 
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Proposition 3.1. Let t[0 : k] be k + 1 distinct real numbers with to = 0 such that 


gaps (t[0 : k]) = g[k] E (0, + 00 )''. 

The distribution of gaps {{X{tQ),X{tk)) has density ^g[k] on where 

k 

^glk]{x[k])= (10) 

TeiS[o,fc] *=i 


In the mono-dimensional case, 


^g{x) = {^g{x) -h $g(-x))l2,>o 


( 11 ) 


and this is also 2$g(x)lx>o when <I>g is even (that is when r = 0 in the stable processes case). 

We may now define a time-homogeneous MC (G^”^[A:] = (G("'\i E [1, kj),n > 0) taking its 
values in (0, -|-oo)^, giving the successive gaps se quen ce starting from an initial one; its Markov 

We will call G^”^ the gaps sequence MC. 


kernel is given by 'k in the sense of Proposition 3.1 
Assume that 


-( 0 ) 


is a r.v. which possesses a density fk on (0,-|-oo)^. The density of g[,^^ is 
OPfc(/fc) where Op^ is the following integral operator (which sends fk onto Opi.{fk)), where for 
any x[k] E 


Opk{fk){x[k]) := j ■ ■ ■ j fk{g[k])'ifgik]ix[k])dgi.. .dgk. (12) 

Of course, if one considers a case for which the iterated process converges in distribution, 
l(-)[k] = (/W(ti),.. -J^^Htk)) ^ I[k] = (/(ti),... ,/(4)) 

n 

then the associated gap MC (G(”)[A:], n > 0) converges too since the map x[0; k] —)• gaps (x[0 : k]) 
is continuous. The converse is false but not that much: the gaps sequence characterises 
the points relative positions. An additional information is needed to recover their posi¬ 
tions: somehow the distribution of the translations which sends gaps {I{ti),i E [[0,A:]]) onto 
{I{ti),i E [0,/;:]]}, and the distribution of the permutation which provides the distribution of 
{I{ti),i E [O, A:]]) knowing {I{ti),i E [0, AJ}. A simple but powerful trick, discussed at several 
places in the paper is the following : we are able to pass from the gaps sequence MC to the 
usual one if instead of {k^\ti),i E [[1, AJ), we study {k'^\ti),i E [O, k}) instead, where to = 0. 

We can sum up in two slogans the relative importance of the iteration of the initial process 
X with respect to the gap MC: the proof of convergence is easier for {k'^\ti),i E [[0,A:]]), but 
the behaviour of G(”)[A:] is easier to understand, and its distribution in the case of Brownian 
processes is tractable. 

3.2 The iterated process evolution 

Any sequence t[0 : k] such that to = 0 can be encoded by the pair C[f\ := [ 9 [A:],t] formed by 
the gaps sequence of t, and the “labelling permutation” r E 5[[0, kj, so that 

h = 5r(i)-5r(o)A e [[0,A:]]. (13) 
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Of course, thanks to (13), the decoding t = C ^{g[k],T) is well defined too (taken to = 0). 
Follows from (13) again, that t,—i(j)is non decreasing in i and then for any i we have 


f'r 9i 9 t{0)‘ 


. , .. .. (14) 

The Markov kernel of the MC n i—>• G [O, fcj) can be made explicit at the level of 

the encodings. Consider t[0 : k] with to = 0 and C[t] := [fl'f/cjjr] its encoding, and t'[0 : k] 
with tg = 0 and C[t'] := [ 5 ^[A:], t'] its encoding. Denote by K the corresponding Markov kernel 
(with transparent convention) which gives the distribution of knowing We 

have 

II, A:l| 


{X{ti) G dt'j^,i G [[1,A:]]} — |x(tj) G dt'^-i^--^,i G 

and then using (j^, (14) and t' = “ 5V'(o) 1°^ * ^ lO) 

Kg[k],T[id9i,---,d9'k),T'] = F{X{ti) £ dt'i,i £ ll,kj) 

= n^Ai; (^C-i(i)) 

i=l 
k 

-.• — 1 


We rewrite more simply as 

k 

Kg[k],r[id9'i: ■ ■ ■, dg'k),T' o r] = 

i=\ 

from what we observe that 

Kg[k\,r[{dg'i, •. •, dg'k),T' ot] = K, 

;HS 

n;" 


(^^9 T'{i) ^ 


(15) 


■g[k],Id\i.d9li ■ ■ ■ idgj^,T^ (16) 

and then the LHS is independent of r. Of 
positive giS, g['s. 


course, all of this is valid for r, r' G 5[[0, fcj, and for 


3.3 Asymptotic independence of labelling permutation and gaps sequence 

We explain now why in the encoding Markov chain (C'[/^”^[/c],re > 1), the gaps sequence 
“becomes progressively” independent from the labelling permutation as stated in the main 
convergence theorems of the paper, where this appears under the form of exchangeability of 
the limiting distribution 7 ^. The asymptotic exchangeability can be proved directly (see [7] 


or the end of Section 5.3). It is somehow quite complex since it relies on the convergence of 


{I^^\ti ),..., I^^\tk)) to a limit independent of the tj’s, and the proof relies on some (classical 
but) involved estimates. 

We present here another argument which makes this more apparent and which we think, 
can be of some interest if ones tries to iterate some processes for which the arguments developed 
in Section 15.31 fail. 
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It is a coupling argument. For a fixed pair {g[k],T), consider (using 


Kg[k],r[^dg[,...,dg'f,),T''] = niiniFg[fc]_^[(dffi,..., r'o r] 

= ^^p^Kg[k]jd[{dg'i,---,dg'k),T'] 


the “minimal flow” going to [{dg'-^,... ,dg'j^,T'] from {g[k],T)^ minimum taken on the t' G 

In general ^ is a defective Markov kernel. Since it does not depend on r", the marginal 
restriction of ^ to the permutation labelling, is the uniform distribution on 5[[0, kj. 
Therefore, Kg^k],T[id9i, • • •) dg'/^), r"] possesses a simpler form: 


Kg[k],TWi,---,dgk),T"] 


^g[k](.dgi, ■ ■ ■ ^dg^) 


lr"e5[0,fc] 

(A: + l)! 


(17) 


where k is a defective Markov kernel on Let 


q{g[k]) = K ( g [ fc ])( M +^) 


be the total mass of and of (Notice that the cases treated in the present paper, 

for g[k] G (0,cx))*^, q{g[k]) > 0.) Now set 


K 


gl],ri^dg[,...,dg'j^),T' ot] 


_^[2] 

^9W\,r 


i^g[k]{dg[,...,d9'k) / jk+iy. 
9{9[k]) 

1-9(9W) 


(18) 


so that is indeed a Markov kernel. It is easily seen that the initial kernel K can be 
represented as 

Kg[klr = q{g[k])K^gl-^^^ + (1 - g(ff[/c]))iL'^i]_^, 

which is the core of our coupling; to sample the MC C'[/(”)] from {g[k\,T)^ hrst, sample a 
Bernoulli random variable with parameter q{g[k]). If it is 1, then use the kernel else the 
kernel If the kernel is used, the new value (G[A: +1], r^+i) has the following property: 
Tfc_|_i is uniform and independent from G[k + 1], which has distribution Kg\^k\{-)/q{g[k]). 

Then as soon as a transition is used the labelling permutation and the gaps sequence 
become independent, and this independence carry on, since by K the labelling permutation 
evolves somehow independently from the current labelling permutation (and it evolves by 
product, see p^). It remains to say some words about the frequency of these renewal events: 
letting = (G[A:]("') , Tn) the successive values of the encoding chain, one sees that each 

time the renewal probability is q{G[k]^'^'^). To get renewal with probability one in the sequence 
C[/W] we need not much: continuity and positivity of the kernel on each compact, and 
tightness of the sequence 


4 Iteration of Brownian processes 

This section is devoted to our results concerning the iterated BM ad libitum, iterated reflected 
BM ad libitum, and nth iterated BM. We will consider iteration of standard linear Brownian 
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motion, but using Remark 5.2 iteration of Brownian motions multiplied by a constant can be 
studied as well. 

We start with a key point relative to the description of the Markov kernel of the gaps 
sequence MC when X is a BM (but many of what follows is valid for more general Gaussian 
processes). In this section, is the density of the centred Gaussian distribution with vari¬ 
ance g. We denote further by Exp[A, x] = Xe~^^lx>o the density of Expo[A], the exponential 
distribution with parameter A. Let MEX^ be the set of probability measures on having a 
density of the form 


f {x[k]) = njExp [Ai,Xi] d^(Ai, • • • , Afc), x[k] £ 


(19) 


\i=l 


where /i is a general probability distribution on called the parameter law of /. In other 
words, the set MEX^ is the set of mixtures of product of exponential distributions. The key 
result in this section, valid only in the Gaussian case, is the following proposition. 

Proposition 4.1. For any k >1, Op^. is linear on MEX^, and then MEX^ is stable by Op;.. 

Proof. We start by the one-dimensional case for which ( [IT| ) holds. 

Let /i(x) = Exp[A, x], and let us find Opi(/i)(x) by computing its Fourier transform 


FTo{a)= / e 
Jx>0 


/ 2^g{x)Xe-^^dgdx. 

'g>0 


This is done in two steps: Opi(/i) is the density of a positive r.v. Z. Hence 


r‘-|-CXD 


FT.ia) = - 


I 2^g{x)Xe ^^dgdx, 
Jg>0 


is the Fourier transform of sZ where e is a uniform random sign, independent of Z. By Fubini, 


one 


finds that it is fg^gXe dg = ^ , which is the Fourier transform of sY where Y 


1 + 


has distribution Expo 




. We deduce from that the identity 


r’+OO 


Exp[A,x]Tg(x)dg'= Exp[\/^,x], x > 0. 


( 20 ) 


In words, Opj^ sends x i—;■ Exp[A,x] on x i—Exp[\/^, x]. 

Remark 4.2. Notice that this implies that Exp[2] is stable by Opj^. This is the result by 
Gurien-Konstantopoulos [7] who proved that R ~ eY where Y ~ Expo[2]. 

Assume k > 1 now. Observe the effect of Op;. on a product of exponential distributions. 
By (20) and (11), one has for any x[k] G (0,-|-oo)*', any r G 5|{0, A:]], the identity 

^ k 1 ^ 

/ ^Yl[^Mci,gi]^g,{^Xx(i))'jdgi...dgk= ^ H l] ' (2^) 

;5lfo.fcll i=i re^lfo.fcll ^ i=i 


Te5[o,fcl 
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An important fact appears here, fact valid only in the Brownian case : one can separate the 
variables Xj’s in the right hand side and let appear a product of independent Expo[4] r.v., 
thanks to the two following identities 


Exp[c, x + x'] = Exp[c, x] Exp[c, x']/c, 

Exp[c,x] Exp[c',x] = Exp[c + c',x]^. 


( 22 ) 


Let us separate the variables, and for this, collect in E^-^i the contribution relative to Exp[.,Xj]. 

Since xi x-^-Q-p then + • • • “f !))• 

Let Er^i = {j : x* G |Ax^q)|} = {j : min(r(j),r(j + 1)) < i < max(r(j),r(j + 1))} be the 
sequence of indices j such that x* appear in |Axt-q)|. Further, let 


WT{c[k]) = 

1 TT \/2cj 

2" FrM 


and Fr{c[k]) = (FV,i(c[fc]),i G [[l,fc]]) where 



FrA4k]) 


(23) 




As a consequence of the previous discussion. 



Lemma 4.3. If f is the map f{x[k]) = 

1 Exp[Ai,Xj] for some fixed A[/c] G (0, -boo)^, 

then 

OPk{f){x[k]) = 

k 

Wr{X[k]) Exp [FV,i(A),Xj]. 


TecS[o,fc] 

4=1 


Of course, this ends the proof of Proposition 4.1 

□ 

In the 2-dimensional case, the 6 functions FV and weights are the following 


^( 0 , 1 , 2 )(c 1 ,C 2 ) = (si, S 2 ) , 

W^( 0 , 1 , 2 )(c 1 ,C 2 ) = 1/4, 

(24) 

•P(0,2,1)(c1,C2) = (si,Sl -b S 2 ) , 

W^(0,2,1)(c1,C2) = 1/4 S 2 /(si-bS 2 ), 

(25) 

■^(l,0,2)(cij C 2 ) = (si -b S2, S 2 ) 5 

^^(1,0,2)(c1,C2) = 1/4 Si/(si-bS 2 ), 

(26) 

-^(1,2,0) (Cl ) C 2 ) = (S2, Si -b S 2 ) 5 

^(1,2,0)(C1,C2) = 1/4 Sl/(S1+S2), 

(27) 

■^(2,0,l)(ci) C 2 ) = (Sl + S2, Si) , 

^i’(2,0,l)(ci,C2) = 1/4 S2/(si-bS2), 

(28) 

^(2,1,0) (Cl, C 2 ) = (S2, Si) , 

^^( 2 , 1 , 0 )(c 1 ,C 2 ) = 1/4, 

(29) 


where for short, we have written Sj instead of\/2ci. We now pass to the consequences in terms 
of iterated BM ad libitum, reflected BM, and in the case of iterated BM. 


4.1 Iteration of BM ad libitum 

Proposition 


1.1 


ensures 


the convergence of G [[1, /cj) to I[k] for any distinct and non 

zeros tj’s, as well as the exchangeability of the limit. Hence gaps(/[fc]) is the limit of the gap 
MC, and the limit of this MC does not depend on the tj’s. The gaps sequence is not sufficient 
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to describe I[k] even up to a permutation (which would be uniform by exchangeability), since 
the gaps sequence determines the set of elements of the sequence up to a translation. We 
present a simple trick which allows one to pass this (apparent) difficulty. 

Consider I[k + 1] a ^ik+i distributed sequence. Take U a r.v. uniform in [I, A: + Ij inde¬ 


pendent from the Jj’s. By Proposition 1.1 


J\k -\- 1] {li — Ijj,i G [[1, k -\- 


(30) 


is a random sequence with one zero entry (with uniform position), and the rest of its entries 
has the same distribution as (Ij, 1 < i < k). Moreover, gaps {J[k -|- 1]) = gaps {I[k + 1]) since 
translations conserve gaps sequence. 

Denote by 7 ^ be the distribution of gaps {I[k + 1]). The following proposition, consequence 
of the previous discussion, allows one to get using 7 ^. 


Proposition 4.4. Consider {Gi,i G [1,^:]]) a random vector distributed according to 7^, U a 
uniform r.v. on [[0, A:]] and r a uniform random permutation taken in 5[[0, A;]], all these r.v. 
being independent. The following identity holds — Gu, 0 < i < A:) = J[A: -|- 1]. 

It remains to describe 7 ^. The case A: = 1 is a consequence of Proposition [m see also 


Remark 4.2 For k > 2, this can be obtained by looking at the limit of the gap MC, starting 


with some initial positive gaps sequence g[k] since, the gap MC inherits from the initial chain 
(the iterated BM) the property to possess a limiting distribution, independent from the starting 
point. Clearly, the initial sequence can be taken random (with values in M'*'^), for example, 
one can start with some independent exponential r.v. with parameters Ai,..., A^... and this 


is what we will do since Proposition 4.1 and Lemma 4.3 allows one to control exactly the 
evolution of the distribution of the gaps sequence MC in this case. 

Finding the limiting distribution in this case amounts to hnding the fixed point of Op^.. 


From Proposition 4.1 and Lemma 4.3 one sees that Op^ sends an element of MEX^ on a 
weighted sums of elements of MEX^, where the total weight is 1: it is a Markov kernel. It 
can be better understood if instead of seeing the action of Op^. at the level of functions, it is 
seen at the level of the parameters (the parameters of the involved exponential distributions): 
consider a (discrete time homogeneous) MC {Z^'^')[k] = ..., Z^'^\k)),n > 0) dehned 

on whose kernel Q is dehned, for any Borelian A of and A [A:] G by 


Q{X[k],A)= ^ w^r[fe](A)(5FRA[fc])(^)- 

TecS[o,fc] 


In other words, 


jp = FV(A[A:]) | Z^^\k] = A[A:]^) = WT{X[k]) for any r G 5[[0,j 


(31) 


(32) 


If Z("')[A:] ~ n, denote by nQ the distribution of Z^'^~^^'l\k\. We can sum up the preceding 
consideration as follows: 


Proposition 4.5. Assume that the gaps sequence {G^'^\i),i G [[1, A:]]) at time n = 0 is a mixture 
of exponential distributions with density go and parameter law then is the parameter 

law o/Op;j(/o). More generally is the parameter law of fn = the density of 

(GW(i),i G [1,A:1). 
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We now conclude by discussing the asymptotic behaviour of the parameter law MC. 

Proposition 4.6. (1) Q is ergodic in [2, 2k‘^]^ (meaning that for any having its support in 
[2,2k‘^]^, converges weakly when n —>■ +oo to a distribution Vk independent from 

(2) The probability density g whose parameter law is is solution to Opi^{g) = g. 


Proof. We prove the two statements. Let M.{S) be the set of probability measures with support 
in S. Since the compact set [ 2 , 2k‘^]^ is stable by any LV, then At([ 2 , 2 A;^]^) is stable by Q. Take 
I/® in A1 ([2, 2A:^]^), and being the parameter law of some function /q. Hence, the sequence 
:= > 0 ) possesses an accumulation point in the compact M-{[2,2k‘^]^). 

(n) in) 

Consider a converging subsequence, still denoted . Recall that v(, is the parameter law 
of fn '■= Op^^\fo). For any fixed x[k], the map X[k] —)■ 0^=1 [-^ i,Xi] is bounded continuous 

on therefore —?■ implies that for any fixed x[k] G ( 0 , + 00 )*^, 

„ k „ k 

fn{x[k]) = j Y\^^p[>^i,Xi\du^^\x[k]) ^ f{x[k]) := I Exp [A*, x,] dr'A:(A[A:]). (33) 

i=l i=l 


The fact that / is a density can be checked by Fubini. Denote by rjn the distribution on 
whose density is fn and by 77 the one whose density is /. By Scheffe’s theorem, the simple 
convergence (33) implies the convergence of pn to r]. 

This implies rj = jk (by uniqueness of the limit of the gaps sequence Markov chain), and 

/ is the density of 7 ^. We must 


then / coincides with lim,i Op^”'^(/o) 


By Proposition 


4.5 


add that a function / in MEX possesses a unique parameter law, which implies that 
possesses a unique accumulation point, and then converges in distribution. The uniqueness of 


the parameter law comes from (19), where one sees that if n is the parameter law of /, then / 
is the Laplace transform of the measure fntiAi)KA[fc]). □ 


Remark 4.7. Take a bounded continuous function / : —)■ M. Our representation of the 
gaps sequence distribution of IBM permits to calculate K{f{G[k])) under 7 ^ and to give a 
representation using Uk only: 


IE(/(G[fc])) = [ f{x[k])d'ykix[k]) 

jRi 


(34) 



duk{X[k]) 


(35) 


hence, it appears clearly that E(/(G[A:])) can be computed thanks to the parameter distribution 
Uk only. More generally, using Proposition |4.4[ one can use this formula to compute E(/(/[/c])) 
to, which can then also be expressed in terms of nj. only. 


MCs with kernel such as Q, that is, which relies on successive applications of a functions 
Ft, where F^- is taken at random in a set of functions F = (TV, r G 5|{0, kj) depending (or not) 
of the current position, are called iterated function system (IFS) in the literature [3l|2l[9]. 
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Figure 1: The support of 1^2 computed by a program. 


Here since 0^*^^ := [2, is stable by all the F-j- (for r G 5[[0, A:]]), it is easily seen that for 


©i”’ - 


U ^-(© 

re-s[o,fcl 


k ) 


the sequence (0^, , n > 0) is a sequence of non increasing compact sets whose (non empty) limit 

(n) 

is a compact Qk- Using the portmanteau theorem and the fact that n 1 —)• 0), Ms decreasing for 
the inclusion partial order (see Figure for a representation of 02) we can establish that for 
any A: > 1, 0^ D Support(i^fc). 

Lemma 4.8. 02 = support(j/ 2 )- 


Proof. For A: = 2, it is easily seen that all the Fj- (given in (23)) are contracting in equipped 
with the Euclidean distance. Following classical theorems (e.g. Hutchinson [121 section 3]) it 
turns out that 0^”^ converges to &k for the Hausdorff metric for any starting set 00 C [2,8]^ 
(and not only from [2,8]^ as stated above). In particular, imagine that 02°^ = {(2,2)}, and 
that the starting measure is = <^( 2 , 2 )- Recall that —>■ 1/2 (since the convergence of 

i/{^) j /2 holds for any starting distribution whose marginals own no atom at 0). 

Take any x G 02 , any e > 0. By Hutchinson’s result, for n large enough r\B{x, e) 7 ^ 0 , 
which means, taken into account the positivity of the that v^^'>Q'^{B{x,£)) > 0: some 
mass has been transported in a neighbourhood of x in n steps from (2, 2). This is a first step 
in our proof that x G Support(02). Now, observe that for any p > 0, there exists m > 1 such 
that 

^oT2([2,8]2)cR((2,2),p), 

implying that m iterations of Foq ^2 (see ([23^), bring back all the mass (that is 1) in a neigh¬ 
bourhood of (2,2). The probability to proceed to these iterations of Foq ^2 is positive (since 
inf(ci C 2 )e[ 2 , 8 ] 2 'W^ 0 ,i, 2 (ci, C 2 ) > 0). Now, since all the functions Fr are uniformly continuous 
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on the compact [ 2 , 8 ], for p small enough, any distribution with support included in 
B{{2,2),p) will also satisfy i'^^^'Q'^{B{x,2e)) >0. □ 


Remark 4.9. Hutchinson m section 3] characterises the set © 2 : it is the closure of the set 
of hxed points of the functions o ■ ■ ■ o Fr^,m > 1, r* G 5|]0, 2]|). 

Using [2] and some analysis, for fe = 3, the IFS with place-dependent probabilities is 
contracting in average for the ||.||2 distance. This can be proved by computing the Jacobian 
matrices Jr{c[k]) = of the Ur’s, and by proving that their norms A^r(c[A;]) := 

supp^o IIP'^t(c[A:])|| 2 /||/o ||2 satisfies r(;T-(c[A:]) log(A^T-(c[A:])) < 0 (this can be proved by taking 
first some bounds on the Wr, and then using the log(A^.r(c[A:])) < log(A^T-(2,..., 2)). From 
Theorem 1.2 in [2], 1/3 is of pure type, atomic, or absolutely continuous. 

We think that the same results can be proved with additional work for k = A, but for k > 5, 
other methods should be involved since the average contraction property seems to fail. We were 
not able to find in the literature any general results allowing one to prove the identification of 
Bfc with Support(/rfc) or to compute the Hausdorff dimension of this support. We conjecture 
that for any k > 1, Support(r'fc) coincides with 0^ , and that the Lebesgue measure of this 
support (or of 0 ^) is 0 . 

Now, we describe the distribution of the gaps sequence of the IBM thanks to I'k- 


Theorem 4.10. Let k be an integer larger than 0. If {Gi,i G [[l,/i;]|) is a random vector with 
distribution 7 ^, then 

{Gi,i€ll,kj) {GiEuG2E2,...,GkEk) (36) 


where the Ei’s are i.i.d., Exp[l] distributed, independent from G[k], a random vector of law Vk- 


According to this theorem and Proposition 4.6, we may deduce the following multivariate 
stochastic order bounds for 7 ^, which somehow, describe the repulsive-attractive property of 
the gaps sequence. 


Proposition 4.11. Let k be an integer larger than 0 and (Gj, i G [[1, A;]]) a random vector with 
distribution 7 ^. For any bounded increasing function h : —?• M 


'E{h{Ei/{2k‘^), I < i < k)) < K{h{Gi, 1 < i < k)) < 'E{h{Ei/2,1 < i < k)) 


where the Ei are i.i.d. random variables Expo[l] distributed. 

We can add here that the bound 2k‘^ is not tight (even in the case A: = 2 as one can see on 
Figure [^. 


4.2 Iteration of reflected BM ad libitum 

In this section A = |R| is the reflected BM (RBM), and = A„ o • • • o Ai the nth iterated 
RBM. 

Proposition 4.12. Let to = ,tk be some non negative distinct real numbers. The 

sequence 

(l(-)(t,),iGlO,A:l) ^(0,/i,...,4) 
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where {Ii,i £ [[1,A:]]) is invariant by permutation and independent from the U’s and takes its 
value in (0,+cx))^. Moreover we have Ii ~ Expo[2]. 

Hence, the gaps sequence MC (gaps i G [[0, fc]]),n > 1) converges and its limit 

{Gi, i G [[1, /cj) := gaps (0, /i,..., If-) determine (0, /i,..., lyf): for a uniform permutation r G 
5[[1,A:]] independent from {Gi,i G [[l,fc]]), 


I[k]^^ {G,^^,i£ 11, 


The proof of this proposition can be adapted from the proof of Theorem |5.1 
results rely on the ergodicity of the Markov chain iyl^'^\ti),i G 


__^ All these 

[[0,A:]]). The estimates needed 
to deal with the reflected Brownian case can be simply adapted from the simple Brownian 
case. 

We now describe the limiting gaps sequence using again the MC at the parameters level. 
First, the Markov kernel of the iterated BM has a density. Set, for any g,x,y > 0, 
Mg(x,y) = F {Bt+g G dy \Bt G dx). We have, by Andre’s reflection principle, 


Mg{x, y) = {^g{y -x) + ^g{y + x)) lg>o 


(37) 


where ^g is the density of Bg. 

The gap MC kernel can be described too adapting consideration of Section (in words, 
0 stay at the left). Starting with some gaps sequence gaps (to = 0, ti,..., tfc) = g[k], we will 
have gaps (^(to) = 0) -^(^i); ■ ■ • = x[k], if, for the same notation as in (§, for i G [I, A:]], 

X{ti) = Xr^j) for some permutation r G 5[[1,A;]] (instead of |10,A;]] for the iterated BM). We 
then have in this case a solid link between the Markov kernel of the gaps sequence and of the 
initial chain, since tj = 9i- We get in this case 


k 

^ 9 [k]ix[k]) = E n 1)) ^'r(i)) lxi>0- (38) 

re-S[l,fcl i=l 


Therefore, modifying a bit 

21 

), 

one can still see that MEX is stable by the MC with kernel 

Tfc. One observes using ( 

20 

, (11 

), 


. k 

/ n 9i]'^g[k]ix[k])dgi...dgk 

R+'= 


+ Exp 


|3^T(i) T ^T(i— 1)1 


2fc 'y y IT Exp |®r(i) ®T(i—1)1 

TecS[l,fc] i=l 

After expanding this product, one can again put together the elements “containing a given” 
Xi- Using the same considerations as those below Remark |4.2[ this is also 

' E E n Exp \x^(^Pi - Exp ly^i, |x.r(j) + a;r(i-i)lj 


2k 


re-S[i,fc] Dc[l,fc] ieD 


ieZD 


This formula is the analogous in the case of RBM to that on the BM, (21). Let E-r^i as defined 
in Section]^ and 


K,i = {j ■ Xi G k^(j)|} = {j : r(j) > i} 

K,i = {j ■ Xi G |T^q-i)|} = {j ■■ T{j - 1) > i}. 
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Set 


and 


Again let 


FT,D,i{c[k\) — V^^Cj 

3&D,j&Er,i jeCDjeE;^ ieCDjes;'^ 


w^r,D(c[A:]) = 


2 ^ -^r,D,i(c[A:]) 

-^r.D = I <i <k). 


Similarly to Lemma 4.3 we have 


Lemma 4.13. If f is the function f{x[k]) = 0*=! Exp[Aj,Xj] for some A[fe] G (0,+oo)^, then 


Opkif){x[k]) = Y1 Y1 'Wr,D{>^[k])YlExp[Fr,D,iW,Xi]. 

Dc[i,fcl Te5[o,fcl 


2=1 


As in the iterated Brownian motion case, to this operator one can associate a Markov chain 
with kernel Q (defined as (31)) at the level of the parameters (see also (|32[)). Again, the 


Markov chain stays eventually confined in a compact region of (the compact [2, 
is conserved by each of the Ft^d)- By the same considerations as that of SectionProposition 
4.6 holds for the present case (with [2, instead of [2,2A:^]^). The analogous of Theorem 

4.10 holds too, for Vk the fixed point of Q, and Proposition 4.11 too, with 18A:^ instead of 2/c^ 
(again 18/c^ is not tight). 


4.3 nth iteration of the BM ad libitum 

In the literature, the standard iterated Brownian motion corresponds to our process It 
has been deeply studied. It permits to construct solutions to partial differential equations m- 
Burdzy studied some of its sample paths properties [5]. Lot of results have been obtained 
around its probabilistic and analytic properties, see [a 12 El El da US] and the references 
therein. The nth IBM permits to construct solutions of differential equations |15j . but they 
are less studied, only [2 mentioned that his result can be extended to nth IBM. As far as we are 
aware of, there are no result concerning some description of the finite dimensional distributions 
of this process. In the sequel, we show that our gaps point of view allows one to give (a non 
trivial) description of them., but sufficiently simple to make some exact computations for small 
values of n and k. 

Let n > 1 be fixed, as well as (to = 0,ti,... ,tk) some distinct numbers. The aim of this 
part is to describe the distribution of (I^”^ {ti),i G [[0, A:]]), where = Bn o - ■ - o Bi, where the 
Bi are i.i.d. two sided BM. 

Start with formula 
4>g is the Gaussian 


We built our reflection on the considerations presented in Section 3.2 


(15) which expresses the encoding Markov chain kernel. Here, of course, 


density. Again, by (20) 




2=1 


g[k],T{{9i, ■ ■ ■, 9k):T' ° T)dgi...dgk = ^ H I 

2=1 


( 39 ) 
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Thanks to (22), we can again rewrite the right hand side to code the evolution on the parameter 


space. Setting m* = min{T-_i, r-} and M* = max{r'_;^, r-} we get = 9i+mi - ^ b'm,- 

Once again, collect the different contribution: set Ej{T,T') = { i : j G \mi + l,Mj]]} the set of 
indices i such that g'^ contributes to The RHS of (39) rewrites 


Wr 


,(A[ti)n Exp j (A), g'jj 
i=i 


where 


WrAm) = ^ n ^ and T.,.^,(A[/^]) = 

- r,r')A\i ieEJr.r') 


r=i Ei 


J = i ^'l■&Ej{T,T') V —-I 

Consider MEX^ the set of measures that are mixtures of distribution on x 5[[0,fc]] of the 
type (nt, expm) X 6r where 5r is a Dirac on a permutation r. The previous considerations 

show that the kernel K operates linearly on MEX^. It sends (OiLi Exp[Ai]^ X 5r on 


'y ^ Wr,T' ) I Exp j ( A[/c] ), Q'j j I X 6-j-'oT' 

t' \i=\ ) 

This can again be written at the parameter level under the form of a time homogeneous MC 
on X 5[[0,/cj, which, starting at time 0 at position (A[A:],r), takes at time 1, the value 
((FT-yj{X[k]),l < j < k),T' ot) with probability WTy{X[k]) (for any r' G 5[[0,A:]]). Denote 
again by Q the corresponding kernel. 

This explicit description allow computations for small values of k and of n. Recall at the 
beginning of Section 3.2 the decoding map C~^ . Finally denoting by the expectation 

when the initial encoding distribution is (nr=i Exp(A.)) X 6r, we find 

Ea[H,.(/(/<”'W)) = E / (40) 

X J^^f{C~^iy[k],T' ot)) dyi...dyk. (41) 


The LHS appears to the Laplace transform of (ti),{tk)) with respect to the initial 
gaps sequence, and then it characterises the distribution. This is not a simple description, but 
we think that it is the simplest representation of the finite dimensional distribution of the 
iterated BM one can find. 


5 Stable processes iterated ad libitum 

5.1 Main results 

In this section, we consider ind^endent two sided-stable processes X\,X 2 ,..., with parameters 
(a,a,r} as defined in Section and there successive iterations = Xn o ■ ■ ■ o Xi. In this 
section, is no more the Gaussian density but the density of -^ 1 ( 15 :). 
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Two sided stable processes possess independent and stationary increments, as well as a 
scaling property which makes their iterations very similar to that of BM (general Levy processes 
seem more difficult to handle because of this lacking scaling property). Here are the convergence 
results we get for iterated stable processes , as described in Section 

Theorem 5.1. Assume a G (0, +oo). Take fe,n > 1 and some non zero ti,... ,tk- Set, 

1. When a < 1 and any r, for any t > 0, does not converge in distribution in M. 

2. When 1 < a < 2 and |r| > 1 then does not converge in distribution in M 

3. When 1 < a < 2 and |r| < 1 the MC converges in distribution. The limit 

distribution fik does not depend on the U’s and is then exchangeable. For I[k] ~ pk, the 
equality I [k] = (X(Ii),--- ,X{Ik)) holds. Moreover, {I 2 — h, ■ ■ ■, Ik — h) ^ Tk-i 
When r = 0, under pi, Ii = eni>o where the X(1)W ’g are i.i.d. copies of 

X(l) and e is an independent uniform random sign. 

Remark 5.2. If X and X' are two stable processes with parameters (a, 1,0) and (a, < 7 , 0 ) for 
some a G (1,2] and <7 > 0, then X' aX. The successive iteration of (a, 1,0) and (a, < 7 , 0 ) 
stable processes and limits (if any) can be compared by a simple coupling, but the property 
fails when r 0 . 


Remark 5.3. Notice that since Pk+i is exchangeable, the rank of Ji in I\k + 1] is uniform. 
Therefore the (random) number of indices ff{j : 2 < j < k + 1, Ij — Ii > 0} is uniform in [[0, kj. 
In other words, if I[k] ~ pk the rank of 0 in the list (0, /i,..., Ik) is uniform. 

Lemma 5.4. ITTien 1 < a < 2 and |r| < 1, the MC {G^'^'i[k],n > 1) converges in distribution, 
and the limit distribution does not depend on the initial non-zero state. 

Proof. Take some gaps sequence g[k] G M*^. They are the gaps sequence of some non zeros 
and distinct times to,ti,... ,tk. Start with /(^^[fc + l] = (ffi, • • • Tfc)- Since /(")[/c + l] converges 
in distribution to a limit independent from the ti (Theorem |5.1[ 3)), then the gaps sequence 
MC [k] := gaps (/(^) [k]) too. □ 

Consider the case k = 1 and X a symmetric stable process for some a G (1,2] and r = 0. 
Assume that the gap at time 0 is distributed as G, at time 1, the gap will be |A(G)|. Then 

from the equality G |A(G)| = G^/"|A( 1 )| we infer 

G^=^n|X(l)«p, (42) 

i>0 


where the A(l)(*) are i.i.d. copies of A(l) (the complete argument can be adapted from Section 


5.3). When r 0 there are not any such simple formula. 
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Let cj) be the density of G as defined in (42). We have Opi((/>) = cf), and (p{g)'ifg{x)dg = 
</>(x) is an identification between the densities of |X(G)| and of |G|. Using that 4>[c,g] := 

is the distribution of cG, we get that Jgy^(j)[c,g]^g{x)dg is the density of \X{cG)\ 
c^/^\X{G)\ c^/'^G, which density is (/i[c^/“,x]. All of this can be summed up in 


/ ff >0 


(t)[c, g\^g{x)dg = 0 [c^/“,x]. 


(43) 


Let v[c\ the distribution whose density is (f)[c ,.]. The map Op^ sends v[c\ onto and is 

a linear application on the set of mixtures of distributions z/[c]. This property which in the 
Brownian case allowed us to prove that Op^ was linear on the set of mixtures of product of 
exponential distributions can not be extended here. The reason is that, to separate the variables 


in (21), we used (22). This important property holds only for exponential distributions, and it 
turns out that in the stable case, product measures of the form 0^=1 are not sent by 

Op^ on mixtures of measures of the same kind. We were not able to find a family of measures 
on which Op^ would operates simply but the discovery of such a family would be an important 
step for the identification of the distribution of I\k]. 


5.2 Occupation measure in the stable case 

As stated in Proposition [T^ Curien and Konstantopoulos [7] obtained some information about 
the occupation measure of the iterated Brownian motion ad libitum. In the stable case, when 
convergence holds, the family of limiting distributions gk are consistent, and since, they corre¬ 
spond to distribution of exchangeable vectors, by Kolmogorov extension theorem together with 
de Finetti representation theorem, there exists a random measure //, so that for any k > 0, gk 
is the distribution of (Ui ,..., Uk) i.i.d. random variables taken under g (this is explained in 
the Brownian case in [7]). 

The main tool used in [7] to characterise the regularity of the density of the occupation 
measure is a paper by Pitt [16] only available in the Gaussian case. We are not able for the 
moment to get a similar result in the stable case, and then we renounce to go on our research 
in this direction. In view of Figure]^ we may expect that for some small parameters a in (1, 2] 
(close to 1 ), the density of the local time should be not positive on the range of its support. 

In the next subsection, we discuss the finiteness of the support of the limiting occupation 
measure. The proof follows the same structure as that of [T] Prop. 7]. Let P be any two-sided 
real process (in our case P = X, or I). The range of P on [a, 6] is defined by 

Rp{a,b) = sup P{t) — inf P{t). (44) 

a<t<b a<t<b 

In the following, set D = RxiO, 1). 

Lemma 5.5. For any |r| < 1 and a G (1,2], /or almost any t ^0, Rj(n){0,t) converges in law 
to a r.v. A which does not depend on t. Moreover, when r = 0, 

OO 

A‘=’nor (45) 

i=0 

where the Di’s are i.i.d. copies of D. 
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Figure 2: Simulation of the local time of the iterated centred stable processes ad libitum, in 
the case a = 1.2,1.5,1.8 and 2. Each of them is made from an histogram made with a sample 
1 < i < 10®) starting from some fixed position. 


Proof. Let An{t) = inf{/(”)(u), 0 < u < t} and Bnit) = sup{/’^”^(u), 0 <u<t}. When r = 0, 


Rj(„+i) (0, t) 


= sup X{v) — inf ^{v) 

A„{t)<v<B„{t) A„it)<v<Bnit) 

= Rx {An{t), Bn{t)) 

= {Bnit) - An{t)r~" D 


By iteration, we get 

n 

i=l 

where Di are i.i.d. copies of D. Since a > 1 and t ^ 0, t°‘ 1 when n -A oo. Now, we have 

to prove the convergence in law of R? * as n —)• +oo. Write 

n—1 n—1 

log I at ' = X] 1-^*1 • 

i=0 2=0 


By the Doob’s inequality [T71 Theorem II. 1.7], for any /3 G M, 


P [H > x] < P 


sup A(t)| > 
0 <t<l 


X 

2 


2 % [X{1)^] 


(47) 
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But, E [X(l)^] < oo if /? < a. For j3 = 1 < a 


[a Mog \Di\ > z ^] < 


Cste 


which is a summable sequence, since a > 1. By Borel-Cantelli’s lemma, fllLi ' converges 
as n —>• + 00 . This ends the proof when r = 0. 

In the general case, write 

Rj(n+i){0,t) = Rx {An{t),Bn{t)) 

< {Bn{t) - An{t)r~'D + r{Bn{t) - An{t)) 


To prove that Rj(„){0,t) converges, we use Theorem 13.0.1 in [T3j. By (48), 


E[i?j{„+i)(0,t)|i?^(„)(0,t)] - Rj(„)i0,t) < (i?^(„)(0,t))" E[D]- {1 - r)Rj(n)iO,t). 


/E \D] \ 

So ifi?,M(0,t) > f 


— r 


= M, then E [i2j(„+i)(0, (0, t)] - Rj(„){0,t) < -1; else 


E [Rj(„+i) (0, t)\Rj(n) (0, t)] — Rj(n) (0, t) < M“ E [D] + 1 = 6, from what we deduce 

E [i?j(n+i)(0,(0,t)] — Rj(„){0,t) < —1 + bl[Qxi]{Rj(n)i0,t)). (48) 

This proves the ergodicity of (0, t); n > 0) by [HI Theorem 13.0.1(iv)]. □ 


By Lemma 5.5 and the arguments of |7l Section 3.2], this proves that (j) has a bounded 
support a.s. 

5.3 Proofs of Theorem 15.ll 


The main technical point (Theorem 5.1 (3)) concerns the convergence of the MC G [[0 ,1 

in the stable case from which we will derive the other convergence theorem of the paper by 
some slight modifications. 

In the proof X{1) stands for the symmetric part of X{1) so that X{t) = rt + |t|^/"X(l). 

1. Assume a < 1, and r G M. One has + Since 1/a > 1, 

it is apparent that |/(")(t)| should become very large. To prove this, we compare with a 
deterministic geometric sequence c” for (1/a) > c > 1. 

Ef|/(^)(t)| >c”||/(”-^)(t)| > inf E(|rx + x^/“A(l)| > c’^) 

For any x > 

P(|rx + x^/“X(l)| > c'^) = 1-: 




—c — rx 

X 


— sx(i)< 


AjOL 
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and since stable distribution possesses continuous density /i at 0 (see Feller m sec. XV(3)]) 
this is 


> 1 - Ch{0) 
= 1-Ch{0) 


xlja 

2 c" 


r;l/« 


1/a 


2 c^ 


for n large enough and some constant C > 0. We deduce 


> c",Vn > 1) > 0. 

When a = 1, for any r, As the distribution of Xi 

is symmetric with respect to 0, {t), I^^~^\t)x['^ 

from which we get {t){r + x["‘ t OILi Taking the logarithm, 

one sees that does not converge in distribution. 


2. The proof we provide here is valid for any a > 0. In the sequel we assume r > 1 (the 

case r < —1 is similar). For a fixed t, + |l(”'“^)(f)|^/“A(l). For r > 1, 

can be compared with a geometric sequence with common ratio s G (1, r). Write 

> inf F{\rx + x^/^X{l)\>s^) 

V / x>s"-^ 

For any x > 

P(|rx + x^/“X(l)| > s") = 1 -P(-s" < rx + a;^/“X(l) < s”) 

> 1-P(rx + x^/“X(l) < 

= 1 -•*«!) 2 ^17^) 

= i-nxd) 

= 1-P(X(1) > (r-s)x^“^/“) 

> 1 - cs(”-T(i-«) 

for n large enough (we have use that P(X( 1 ) > v) < c'v~^ for some c' and v > s, and that 
r — s is a constant, and the symmetry of the distribution of A(l)). We deduce from that that 

P(|/(’^/t)| >s'‘,Vn> 1 ) > 0 . 

3. The proof of the convergence of we propose is adapted from Curien & Konstan- 

topoulos [7]. 
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The sequence (/^”'^[/c],n > 1) is a MC, and its Markov kernel is given by 

P(y[fc];^)=P((X(yi),...,X(yfc)) E A), 

for any y[k] E any Borelian ^ E As in [7], the Markov chain is aperiodic, and 

irreducible with respect to the p-dimensional Lebesgue measure on We prove that it is 
Harris recurrent (and then possesses a unique invariant distribution), following the elements 
that can be found in Section 5.5.1. Meyn &: Tweedie M- Set for any M > 0 

Sm = E M^, M~^ < \xi\ < M, \xi — Xj\ > M~^}. 

Denote by the density of (A(xi),... ,X{xk))-, and let 

FM{z[k])= min f^[k]{z[k]). 
x[k]&SM 

It is easily seen that Fm is the density of a cr-finite measure ^lm on with total mass 
Cm = /fjfe FM{z[k])dzi ■ ■ ■ dzk > 0 and satisfy FM{z{k\) > 0 for any zi,... ,Zk- This provides 
the following bounds on the Markov kernel of our MC: 

P((X(xi), .. .,X{xk)) E A) > for all x[k] E Sm- 

This is the minoration condition (5.2) in [T3]: the set Sm is /iM-petite. To prove the Harris 
recurrence of the MC, it suffices to prove that for some M > 0, the expected hitting time of 
Sm by lA')[k] starting from x[k], is bounded above for x[k] E Sm- Consider, for x[k] E 

V{x[k]) = U{x[k]) + G{x[k]) 

with U{x[k]) = max{|xj|,z E [[l,fc]]}, G{x[k]) = Ylo<i<j<k 1^* “ Xj\~^A (where xq = 0). The 
potential function V is unbounded on and its drift is defined by 

DV{x[k]) := PV{x[k]) - V{x[k]) = E(H(X(xi),... ,X(xfc))) - V{x[k]), for x[k] E 

We just have to prove that 

XV{x[k]) <-a + hlsM^x), x[A:] E M^. (49) 

We have for any A > 0, 


PU{x[k]) = E ^ max^ 1X3;^^ = E[max |Xa,J + |r| |xi|] < |r|C/(a:) + A^/“E[max |X3,J;)^|] 

and then taking X = U{x[k]), we get 

PU{x[k]) < |r|[/(x[A:]) + C/(x[A:])^/“Ci 
where Gi = E[max_i<s<i |Xs|]. Now, 

PG{x[k]) = 

0<i<j<k 

= X] ® \X^i-xj +r{xi-Xj)\~^/°' 

0<i<j<k 
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We decompose each term in the sum using 
1 


\Xxi-xj + r{xi - 


^Sign(Xj:._j;^.)=Sign(r) ^Sign(X3:)/Sign(r) 

\Xxi-xj + r(xi - Xj)\^/^ 

1 


^Sign(Xj; )=Sign(r) 

< -- _ ^ -- + 


Sign(X^._x)7^Sign(r) 


\Xxi-xj\^/'^ \\Xxi-Xj \ - \r{Xi -Xj)||Vo 

By symmetry and unimodality of the density of centred stable distributions, one has 


E 


Xx; —.7 


+ r(xi < 2E 


Hence 


PG{x[k]) = 2 ^ (|a:i-x,f/“)-^/“E(|Xi|-i/“) 

0<i<j<k 

this last inequality come from which can be viewed 

as an application of Jensen inequality: take W uniform in [[l,m]], f{x) = Since 

/ is concave E{f{\yw\~^^°‘)) < f{E{\yw\~^^'^)) which is equivalent to Zl™ i < 

(m Z2=i ■ We get, by convexity, for some constant Ck and C^, 

PV{x[k]) = PU{x[k]) + PG{x[k]) < Gk{U{x[k]f/^ + G{x[k]f/^) + \r\Uix[k]) 

< G'kV{x[k])^/'^ + \r\V{x[k]) 


which implies 


AVix[k]) < C'kV{x[k])^/° - {l-\r\)Vix[k]). 


(50) 


If x[/c] ^ Sm then there exists i such that \xi\ > M or {i, j) such that \xi—Xj\ < 1/M. In the first 
case V{x[k\) > M and in the second one, V{x[k]) > For M > 1, we thus have V{x[k\) > 

MV“ for X i Sm- The RHS of Q rewrites H(x[A:])i/“(C^ - (1 - \r\)V{x[k]Y-^l^). For M 
chosen such that (I — > max{l,2C'^}, V{x[k]Y^{G/ — (1 — |r|)I/(x)^“^/") < 

For x[k] e S'm, 0 < V{x[k]) < M + {k + < +00 and then AV{x[k]) is bounded on 

Sm, this allows one to prove that ( [i^ holds for G = Sm and M large enough. 

To end the proof, we need to prove the exchangeability of I[k]. The argument is general, 
and present in [7j. Take a G 5|{1, kj and ti,... ,tk distinct and non zeros. By the proof above, 
both G [[l,fc]]) and (/W(t^(i)),i G [1,/c]]) converge to I[k]. So, i G ll,kj) 

converges to I[k] and to (/o-(i))* £ Hence, I[k] G [I, A:]]) for any a. □ 


6 Conclusion 

In the paper, we have presented some results and some tools allowing to study iterated inde¬ 
pendent processes. Our tools are really useful only for processes with increments independent 
and stationary. Hence, the global frame is that of Levy processes. But what we did for station¬ 
ary process could probably done for continuous MC, homogeneous or not. For example it is 
likely that one can get some results on iterated Ornstein-Uhlenbeck processes whose increment 
are simple enough to be controlled. 
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